Dephasing effects on stimulated Raman adiabatic passage in tripod configurations 
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We present an analytic description of the effects of dephasing processes on stimulated Raman 
adiabatic passage in a tripod quantum system. To this end, we develop an effective two-level model. 
Our analysis makes use of the adiabatic approximation in the weak dephasing regime. An effective 
master equation for a two-level system formed by two dark states is derived, where analytic solutions 
are obtained by utilizing the Demkov-Kunike model. From these, it is found that the fidelity for the 
final coherent superposition state decreases exponentially for increasing dephasing rates. Depending 
on the pulse ordering and for adiabatic evolution the pulse delay can have an inverse effect. 

PACS numbers: 32.80.Qk, 33.80.Be, 42.50.Dv, 42.50.Le 



I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) dHi 
is a powerful and robust technique for achieving complete 
population transfer in three-state quantum systems. By 
using two pulsed laser fields population is adiabatically 
transferred from an initially populated state fa , to a tar- 
get state fa via an intermediate state fa. A unique fea- 
ture of this technique is that the intermediate state is 
never populated. This is due to the fact that the sys- 
tem at all times adiabatically follows a dark state, hence 
population losses due to spontaneous emission are sup- 
pressed. 

Apart from being used for population transfer, STI- 
RAP can also be used to create coherent superpositions 
of states in quantum systems in tripod configurations 
@, HL § • The main idea behind this method is the same 
as with STIRAP in A-configurations, with the exception 
that the system now adiabatically follows a superposition 
of two dark states. The interference of these two states 
results in a coherent superposition between two or three 
of the ground states of the tripod system. The exact 
form for this final state is defined by the geometric phase 
acquired by the dark states. 

In addition to the creation of coherent superpositions 
j STIRAP in tripod configurations can be further 
exploited to implement quantum gates 0-0] . Further- 
more, adiabatic passage in tripod systems can be used 
to engineer non-Abelian gauge potentials for ultracold 
atoms The formation of such potentials is made 

possible because during the adiabatic evolution, the sys- 
tem can acquire non-Abelian phases 



11] 



Since the intermediate state is not populated in the 
adiabatic limit, spontaneous emission from this state is 
expected to have no effect on the fidelity. On the other 
hand, the creation of a coherent superposition relies on 
the formation of coherent dark states. Thus maintaining 
coherence is vital for achieving higher fidelities. However 
phase relaxation effects induced for example by elastic 
collisions or laser phase fluctuations can have an adverse 
effect on the fidelity. 



Previous studies on STIRAP in the presence of dephas- 
ing [H|,[l3|, have shown that decoherence can lead to pop- 
ulation losses from the dark state, resulting in a transfer 
efficiency reduction. On the other hand, increasing the 
relative delay between the two laser pulses increases the 
transfer efficiency. This is due to the inverse dependence 
of the transition time with respect to the delay. In a re- 
cent paper by M0ller, Mads en and M0lmer [8|, dephas- 
ing in tripod systems and the effect it has on single qubit 
gates was considered. Using the Monte Carlo wavefunc- 
tion method, they were able to show that the system ac- 
quires complex geometric phases, implying losses which 
reduce the gate fidelity. 

In the present paper, we extend the method used in 
Ref. [IH to study dephasing effects on STIRAP in tripod 
configurations. The method makes use of the adiabatic 
approximation in the weak dephasing regime, where we 
derive an effective two-level master equation for the dark 
states. Analytic solutions are obtained when the Stokes 
and control pulses overlap, whereas for other pulse order- 
ings we make use of numerical simulations. Depending 
on the pulse ordering, similar with or different features 
from STIRAP in a A-configuration are observed. 

The paper is organized as follows. In Sec. |TI]we provide 
a brief introduction on STIRAP in tripod configurations 
and the master equation is introduced. In Sec. IIIII we 
present the effective two-level model for the dark states, 
and derive analytic solutions in Sec. [IV] In Sec. IVl we 
present results from numerical simulations. A summary 
of the results is given in Sec. [Vj 



II. THE TRIPOD CONFIGURATION 

The tripod system is shown in Fig. [TJ The three 
ground states fa, fa and fa are resonantly coupled to 
the intermediate level fa via a pump Q p (t), a Stokes 
f2 s (f) and a control pulse f2 c (£) respectively. The system 
is initially prepared in state fa, i.e. oo) = fa. Such a 
configuration was part of a proposed scheme for creating 
and phase probing coherent superpositions (J- This was 
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FIG. 1. (Color online) The tripod configuration and the three 
laser pulses: fi p (f) (pump), Q c (t) (control) and fi s (t) (Stokes). 
The relative dephasing for each pair of states is depicted by 
wavy lines. 



demonstrated in an experiment by Theuer et al. Q . 

The scheme we examine here uses an adiabatic passage 
method which is an extension of STIRAP in A-systems. 
The population initially placed in state ■01 is either trans- 
ferred to the other two ground states ip3 and ip4 or is split 
between the states tpi an d ip3- The final superposition 
between 03 and ^4, or between 03 ip±, is determined by 
the relative distance of the pulses. The conditions for the 
robust creation of such superpositions is that the system 
evolves adiabatically @, 



where the time dependent mixing angles <j>(t) and 8(t) 
read 



tan(0(i)) 



n e (t) 



and 



tan(0(t)) = 



n P (t) 



(3a) 



(3b) 



We note here, that since the two dark states are eigen- 
states of H{t) with zero eigenvalue, i.e. iJ(t)$i j2 (i) = 0, 
any linear superposition of these two states is also a dark 
state 0, HI- In addition to the two dark states, there are 
also two adiabatic states with non-zero time-dependent 
eigenenergies 



$ 3 (i) =7^1 sin6» + -j=tjj2 
H — j=ipA cos 9 sin 0, 
$ 4 (0 = 7/|^i s[li9 ~ y|^2 



1 



V2 



0>3 cos COS 



(4a) 



V2 



■03 cos 9 COS 



1 



-.tpi cos sin 0, 



with eigenenergies 

*s(*) 



-e 4 (t) = ^n(t), 



(4b) 



(5) 



A. Adiabatic states: Creating coherent 
superpositions 



The interaction Hamiltonian in the rotating wave ap- 
proximation reads 



/ fi p (i) 

n p (t) n a (t) n c (t) , 

n»(t) 

V n c (t) / 



where fi(t) = JO^i) + n 2 s (t) + fl 2 c (t) is the rms Rabi 
frequency. 

In the adiabatic limit, the time-dependent eigenstatcs 
are weakly coupled and this is also valid for the two dark 
states &i(t) and &2{t). Although they form a pair of 
degenerate states, ei(t) = e 2 (t) = 0, the corresponding 
diabatic coupling is always zero, 



(*i(f)|$ 2 (t))=0, 



(6) 



where we take all transitions to be resonant with the 
respective laser field. Furthermore, we assume that 
these are the only allowed dipole transitions. This time- 
dependent Hamiltonian has two dark states $i(t) and 
*a(t) 1 

$i(i) =—=-0i cos# p03(sin#cos0 + isin0) 

v 2 v 2 

p04(sin#sin0 — icos0), (2a) 

v2 

<£>2(i) = ^^0i cos6> =03(sin0cos0 — isin0) 

V2 V2 

p04(sin#sin0 + icos0), (2b) 

v2 



(*i(f)|*i(t)) = -<* 2 (*)|*2(t)} = i^sinfl. (7) 

Because of this time-dependent energy shift, both dark 
states acquire a geometric phase dj [14 1 



i?i = -tf 2 = i9 9 = / di0sin0. (8) 



Thus when the system starts in state $i(t) or $ 2 (i), 
it will adiabatically follow this state and remain in this 
state at all times acquiring a net phase ±$ g respectively 

$l(-oo) ->• e _i ^$i(oo), $ 2 (-00) -> e i,5 3$2(oo). (9) 

As said in the begining of this section, the system is 
initially prepared in state 0i . Then in the adiabatic limit 
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and for 8{— oo) = 0, the system state is a symmetric 
superposition of the two dark states 



The dissipator matrix D describes dephasing effects 



*(-oo) = tpi = -)= ($i(-oo) + $ 2 (-oo)) 
v2 



(10) 



and for t — > oo it reads 



\&(oo) 



1 

7t 



(e- i,5 »$ 1 (oo) + 6^*2(00)). (11) 



The final superposition state in terms of the bare states 
depends on the ordering of the three pulses. Of the many 
possible pulse orderings four different pulse sequences are 
particularly interesting for each produces a different co- 
herent state 0. These pulse orderings arc as follows. 

• The pulses are ordered so that the Stokes pulse 
starts first and ends after the pump pulse, while 
the control pulse is delayed with respect to both 
of them. For this pulse sequence we have the fol- 
lowing asymptotic relations: 8(— 00) = 6(00) = 0, 
<p(—oo) = and <j>(oo) = ir/2. Then the final state 
reads 



^(00) = i/ji cos -dg — 03 sin -d g 



(12) 



• A different pulse ordering is arranged so that the 
Stokes pulse comes first, followed by the control 
pulse and then the pump pulse. For this case the 
asymptotic relations are 9(— 00) = 0, 6(00) = ir/2, 
<f>(—oo) = and 0(oo) — ir/2. Then the final state 
is 



\I>(oo) = —tpssm'&g — 04 cos$ s 



(13) 



• Alternatively one can reverse the order of the 
Stokes and control pulses, i.e. the latter pulse pre- 
cedes the Stokes pulse. Then for t —> ±00 we 
have 6 (—00) = 0, 8(00) = tt/2, </>(— 00) = tt/2 and 
4>(oo) — 0. The final coherent state reads 



^(00) = —■03 cos $ g + 04 sin $ s 



(14) 



• Finally, the Stokes and control pulses can coincide 
in time and precede the pump pulse. Then the 
following asymptotic relations apply: #(—00) = 0, 
8(oo) = tt/2 and <j>(— 00) = 4>(oo) = 7r/4. Because 
<j){t) is constant, we have that (f>(t) = and thus 
the geometric phase is zero. Hence the final state 
is 



*(oo) 



1 

V2 



(03 + 04 



(15) 



B. Dephasing 

In order to model the effect of dephasing in the tripod 
system, we make use of the master (Liouville) equation 



ihp= [H(t),p]+D. 



(16) 



d = -in 



712/012 713/013 714/014 

721P21 723/023 724P24 

731P31 732/032 734/034 

741/041 742/042 743/043 



(17) 



where 7^ = jji are the constant relaxation rates and p 
is the density matrix in the bare basis m , i.e. p mn = 
(0 m |/5|0 ra ). The initial conditions are pn{— 00) = 1 and 

Pmn = f° r mn 7^ 11- 

The derivation of master equations such as the one in 
Eq. (fil)]) is based on the use of the Born-Markov approx- 
imation [l5| . This imposes restrictions on the spectral 
properties of the heat bath with which a quantum system 
interacts [l6| . while the relative coupling strength must 
be weak. Thus, collisions between atoms or molecules 
must be weak [17]] , whereas the fluctuating laser phase 
must be well approximated by a Markovian process [HI . 

In the following section we derive approximate solu- 
tions in the weak dephasing limit and for adiabatic evo- 
lution. Because of this latter assumption we will be using 
the density matrix in the adiabatic basis, defined from 
the following transformation 



P 



RT l pR 



(18) 



where p a is the density matrix in the adiabatic basis 
$j(t), i.e. pf nn = ($ m (t)|p| $„(*)). The rotation matrix 
R is formed by using the adiabatic states as its columns 



(19) 



where R 1 — i?0 Then the master equation in the adia- 
batic basis is 

ihp a = [H a (t), p a ] - ihlR^R, p a ] + R^DR, (20) 

where H a (t) is the Hamiltonian in the adiabatic basis, 

H Q (t)=diag(0,0,e 3 (t),e 4 (t)). (21) 

The second term on the right hand side of Eq. (|20p . 
i.e. R~ 1 R, describes non-adiabatic interactions (off diag- 
onal terms), whereas the two diagonal terms (R~ 1 R)n = 
— (R~ 1 R)22 are responsible for the geometric phases ac- 
quired by the two dark states. 

III. THE TWO-LEVEL APPROXIMATION 

As mentioned earlier, the creation of coherent states 
is the result of interference effects between the two dark 
states ©. This suggests that it is possible to solve the 
master equation (Til)]) approximately by deriving an effec- 
tive two-level master equation. To this end, we make two 
approximations: the adiabatic and the weak dephasing 
approximations. These were used before when studying 
the effect of dephasing on STIRAP in A-systems [13l ]. 
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A. The adiabatic approximation 

The first of the two approximations that we are using 
is that of adiabatic evolution. In order for this to be valid 
we need large pulse areas @, H| so that 



\0(t)\ « \n(t)\, \m « 



(22) 



When these two conditions are satisfied, the system adi- 
abatically follows a superposition of the two dark states, 
whereas the adiabatic states $3(t) and $4(t) are not pop- 
ulated. Because of this the coherences pf- related to these 
two states are negligible. 



B. The weak dephasing approximation 

The second approximation is that of weak dephasing, 
i.e. the relaxation rates 7y are much smaller than the 
rms Rabi frequency Q(t) 

7<i<|fi(*)|, (i,3 = 1,2,3, (23) 



Assuming that the adiabatic approximation Eq. (|22[) is 
valid, then for 7^ > |f2(i)| we would have that jijT 1, 
with T a characteristic time length for the pulse dura- 
tions. This latter inequality corresponds to strong de- 
phasing which would lead to complete incoherent dynam- 
ics that are governed by rate equations [l9|. This justifies 
the use of the weak dephasing approximation Eq. (|23[) . 

Using both approximations we can now simplify the 
analysis by neglecting all coherences that include the two 
adiabatic states <l>3(i) and 3>4(t), 



P, 



0. 



(ij ^ 12,21) 



(24) 



With this the density matrix in the adiabatic basis p" 
acquires the following approximate form 




Pti Ph 

Plx Ph 

p% 3 

fa 



(25) 



Using Eqs. flU]), (HH), ^ and ^ (sec appendix 
we can first show that the population inversion for the 
bright states, w^ 4 (t) = p% 3 — p% 4 , and for the dark states, 
w i2(t) = P11 ~ P221 are both zero at all times, i.e. 



P11 = P221 



P33 = Paa- 



(26) 



Hence, the populations for the adiabatic states can be 
parametrized in terms of a single function s(t), i.e. 



P11 = P22 



P33 — Pa 



1 - 2s(t) 



1 + 2s(t) 



(27a) 



(27b) 



From this we can derive a set of coupled differential equa- 
tions for s(t) and the coherences for the dark states, 
u(t) = ^/2Rc{pl 2 } and v(t) = V2Im{pf 2 }, 

s(t) = - T s (t)s(t) + V2~n su (t)u(t) + y/2Sl m (t)v(t), 

(28a) 

ii(t) =-T u (t)u(t) + (2(j)sm6 + tt uv {t)j v(t) 

+ V2Cl su (t)s(t), (28b) 

v(t) =-T v (t)v(t) + (-2<])sm6 + n uv (t)) u(t) 

+ V2~n sv {t)s(t). (28c) 

The initial conditions are s(t) = —1/2, u(t) = l/\/2 and 
v (t) = 0. The derivation of these equations and the ex- 
pressions for the Rabi frequencies n su (t) , fl sv (t) , fl uv (t) , 
and those for the relaxation rates Tj(t), j = s,u,v, are 
provided in appendix [A] We note here that the effective 
relaxation rates Tj(t) and the Rabi frequencies fiy(i), 
depend only on the relaxation rates 713, 714 and 734, see 
Eqs. (jA8j) and K9]) . 



Fidelity for the target coherent superposition 
state 



In the analysis that follows we will be using the fidelity 
F(t) for the final coherent state *(oo) [20I ] 



F 2 (t) = |(*(oo)|^)|*M)| 5 



(29) 



where *(oo) is one of the four states (JT2J), (HU), JT4]) or 
(|15|) . and p(t) is the density matrix in the bare basis. For 
t — > 00 the fidelity can be expressed in terms of the dark 
state populations and coherences. Using Eq. (fTTj) for 
the coherent state ^(00), and Eq. (|25|) for the density 
matrix in the adiabatic basis and for weak dephasing, the 
fidelity reads 



F 2 (oo) =^(oo) + cos(2tf ff )ReK 2 (oo)} 



sin(2^)Im{^ 2 (oo)}, 



(30) 



where -d g is given by Eq. ([S]). 

To this end, and before proceeding with the solution 
of Eqs. (|28|) . we introduce one more tool that we will be 
using in the following sections. This is the transition time 
T tr (e), and is defined as the time needed for the fidelity 
to rise from a small value F 2 (t e ) = e to F 2 (ti_ £ ) = 1 — e, 
i.e. 



t,. 



(31) 



The importance of the transition time was discussed on 
a previous work on dephasing effects on STIRAP in A- 
systems [l3| . For STIRAP the pulse delay has an inverse 
effect on the efficiency, and this is due to the fact that 
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the transition time is inversely proportional to the delay 
time. 

We should note that the above definition for the tran- 
sition time, does not apply when the Stokes pulse ends 
after the pump pulse, see Eq. (|12p . While for the other 
three possible pulse orderings the fidelity is initially zero, 
i.e. F 2 (— oo) = 0, for the former pulse ordering the fi- 
delity is F 2 (— oo) = cos 2 (?9 g ). For this case, instead of 
using the above definition for the transition time, we will 
be using the following one 

T tr {e) = t!_ e - t 1+e , (32) 

where the time ii+ e is such that 

F 2 (t 1+e ) = (l + e)f 2 (-oo), (33) 

and t\- e remains the same. 

IV. ANALYTIC SOLUTIONS 

A. Equal relaxation rates jij — 7 and Q s (t) = Q c (t) 

A special case is that when the Stokes and control pulse 
coincide, i.e. f2 c (t) = fi s (t). Then at all times we have 
<j>(t) = 7r/4 and (j>(t) = 0. When taking all relaxation 
rates equal i.e. jij = 7 we have that 



n sv (t) = n uv (t) = 0, 



(34) 



and the equations for s(t) and u(t) decouple from that 
for v(t), i.e. 



ci(t) = V2Q su (t)c 2 (t), 

62(f) = A su (t)c 2 (t) + V2n su (t) Cl {t), 

and for v(t) we have 

v(t) = -r v (t)v(t). 

The new variables c\{t) and c 2 (i) are 



(35a) 
(35b) 



(35c) 



ci(t) = s(t)cxp(J r s (i')di' 
c 2 (t) = u(t) cxp ( [ T s (t')dt' 



(36) 



This parametrization is used in order to emphasize the 
analogy to a two-level system. The effective relaxation 
rates and coupling in terms of the mixing angle 9(t) read 



r„(i) = 7 cos 2 (0), 

7 ^3 sin 2 (2(9) 
4 



^SU (^) 

A s „(t) = 



4 

7 ^3sin 2 (2£ 
4 



cos 2 



cos 



7 sin 



(37a) 
(37b) 

(37c) 
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FIG. 2. (Color online) (a) The coupling Q su (t) and the de- 
tuning A su (t) Eq. (J3S]) for Gaussian pulses Eq. (gU). (b) The 
coupling £(t) and the detuning A(t) Eq. (|45|) in the adiabatic 
basis igSJ. The delay is r = T and 7 = T" 1 . 



In Fig. Ufa), we plot the coupling Cl su (t) and the detun- 
ing A su (t) for Gaussian pulses Eq. (|4"4")l . To this end, we 
note that since v(— 00) = equation (|35c[) has the trivial 
solution v(t) = 0. 

In order to connect Eqs. (|35a[) and (|35b[) with those for 
a driven two-level system and exploit this to solve them, 
we make use of the following time-dependent rotation 



where the angle £(t) is 



- cos £ sin £ 
sin£ cos£ 



1 / 2ygn^ \ 

With this equations Q35a[) and (|35b|) take the form 



c_(t) 
c+(<) 



e _(t) £(t) 



c+(t) 



(38) 



(39) 



(40) 
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where the "energies" 



1. Adiabatic following 



e±(f) = ^#(l±sec(20), 



are the eigenvalues of the matrix 



W(t) = 



V2Cl su (t) 



V2n su {t) A su (t) 

The corresponding eigenstates ip±(t) are 

ip+(t) = [sin£,cos£] T , 
ip-(t) = [— cos£, sin£] T . 



(41) 



(42) 



(43) 



Multiplying now both sides of Eq. ([4T)|) with the imagi- 
nary unit, we see that it takes the form of a Schrodinger 
equation for a two-level system. The two states and 
-0+, are driven by a laser field with Rabi frequency £(i), 
while decaying with rates e-(t) and e+(t) respectively. 
The initial conditions for Eq. (|40p can be derived from 
Eqs. (|3H|), ([M]), (JSTcJ) and (j37b|) . and are c+(-oo) = 1 
and c_(— oo) = 0. 

In Fig. [2jb) we plot the "chirped" detuning A(t) = 
e+(t) — e-(t) and the coupling £(i) for Gaussian pulses 
of the form 



(44) 



The "chirped" detuning A(t) and the coupling £(i) are 



(^) 



• t /4tr 



(45) 



where 




(46) 



1 + 5 cosh(a;) — 4 sinh(x) 



From Fig. [2 we see that both £(£) and A(t) resemblance 
the laser field and the frequency chirp for the Demkov- 
Kunike (DK) model [2l|. This feature is exploited next 
to derive analytic solutions for Eqs. (f4*tj)) . To this end, it 
should be pointed out that the time-dependent detuning 
in the original DK model corresponds to real frequency 
chirp, whereas in the two-level system of Eq. (|4U|) it is 
an imaginary chirp, i.e. a time-dependent decay rate. 



The method we use to solve Eq. (|40)) . is the one used in 
Rcfs. [22|,[23| to derive the propagator in three-state sys- 
tems with pairwise crossings. Starting with the system 
initially in state ^>+(oo), we assume that it evolves adi- 
abatically until reaching the crossing point t = 0, where 
A(0) = 0. At this point, diabatic transitions will oc- 
cur and the new system state will be a mixture of the 
two adiabatic states ip±(t). The effect of the crossing is 
expressed in terms of a transition matrix 



U c (0) = 



U—(0) 01+ (0) 
U+-(0) U ++ (0) 



(47) 



where Ujj(0) and Uij(0) (j ^ i) are the survival and 
transition probability amplitudes, for a two-level system 
driven by a laser pulse £ e / / (t) « £ (t) , in the presence of 
a time-dependent spontaneous emission A e ff(t) ps A(t). 

For times t > 0, i.e. after the crossing the system will 
evolve adiabatically. With this the final system propaga- 
tor /7(oo,— oo) reads 

[/(oo, -oo) = U a (oo, 0)U c (0)U a (0, -oo), (48) 

where the adiabatic propagator U a (tf,ti) is 

/ e ftf e-{t)dt q \ 





U a {tf,U) 



e /t/ £+{t)dt 



(49) 



Upon using the above equations with the initial condition 
tp(— oo) = ip+(— oo), the final two-level state takes the 
form 



V>(oo) =C/_ + (0)e / « O ° e - ( * )d * +/ -- e + ( * )d V-(oo) 



(50) 



^+(oo). 



2. Diabatic transitions 

The survival probability amplitude J7++ (0) and the 

transition amplitude U \-(0), are solutions of the 

Schrodinger equation 



(51) 



As already pointed out earlier and shown in Fig. HJb), 
this two-level system resemblance the Dcmkov-Kunikc 
model [2l|, with the only difference that the real fre- 
quency chirp for the latter model is replaced by an imag- 
inary one, i.e. population loss. 

To proceed with the solution of Eq. (foTj) with the ini- 
tial condition c+(— oo) = 1 and c_(— oo) = 0, we utilize 
a hyperbolic sech function to approximate 



&=//(*) = Asech((t - t max )/T eff ) , 



(52a) 
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and a hyperbolic tanh function with a constant term for 
A(t) 



A eff (t) = D + Btanh ((t - t max )/T eff ) . 



(52b) 



where r(x) is the gamma function [26|. The parameters 
a, (3 and S are 



a = AT. 



eff 



— arctan(2-\/2), 

27T 



(58a) 



The time t max refers to the maximum of £(i) Eq. (|45 
and is 



T 2 / 4 

'max = — arctanh I - 



(53) 



The amplitude A for £ e f f (t) is derived from the condi- 
tion i{t m ax) = Lff{tmax), and IS 



.4 



V2T 2 



(54) 



The effective pulse duration T e f f is obtained by requiring 
that both and £ e / } (t) have the same pulse area. With 
this T e ff reads 



T. 



eff 



-arctan (2V2). 



Finally, the condition for obtaining D and B, is that 
for t « t m ax, A(t) A e ff(t). After performing a Taylor 
series expansion for both A(t) and A e ff(t) at the vicinity 
of tmax and keeping only terms of first order in t, we have 



D 



4V67 
^5~' 



and 



(56a) 



(56b) 



Solutions for Eq. (|51[) and for the Rabi frequency 
ieff{t) Eq. ||52g| and for A eff (t) Eq. (|52b| . can be 
expressed in terms of hypergeometric functions, see for 
example Refs. [24|, [H|. For t — ¥ 00 the probability am- 
plitudes are 



B = - 64v ^ 7 arctan(2\/2) ■ 

1257T 



(7++(0) =- 



, r(i + s + /?) 
' r (i + s - s/WTa*) 



(57a) 



/3 



BT, 



eff 



DT, 



eff 



16V^arctan 2 (2V2) 7T 2 
125tt 2 ~ 

2\/3arctan(2v / 2) 7T 2 
25^ ~' 



(58b) 



(58c) 



3. Populations and coherences for the dark states 



After substituting Eqs. (|57aj) and (|5Tb[) into Eq. (|50]), 
and using the inverse rotation 7?. _1 (oo) along with Eqs. 
(|36p . we obtain s(oo) and m(oo). From these we derive 
the real part for the coherence p" 2 (oo) of the two-dark 
states ©, 



(55) ReK 2 (co)} 



.(O)exp 



CulT 2 
At 



(59a) 



and their populations 



p a n (oc) = faioo) = \ + ^[/_ + (0)exp 

(59b) 

where c s = 2.42 and c M = 0.68. The exponential terms in 
the above expressions are contributions acquired during 
the adiabatic evolution, see Eq. ((50)) . Their dependence 
on the different parameters, i.e. the factor 7T 2 jr can be 
easily derived, whereas the two factors c u and c s can be 
calculated numerically, see appendix iBl 

We should recall here that the imaginary part of 
Pi 2 (oo) is zero, because v(t) = 0. Furthermore we note 
that in the long time limit the real part is also zero. Then, 
the fidelity for a coherent state (|30[) . will be proportional 
to the population of the dark states. As we can see from 
Eq. (|59b[) this will decrease for increasing relaxation rate 
7, or will increase for increasing delay r. This is because 
of the inverse dependence of the transition time with re- 
spect to t. The transition time T tr (e) for e = 0.1 can be 
derived from the following expression for the fidelity 



F\t) 



exp 



2 + exp 



'4&T> 

.T 7 ' 1 



(60) 



and is 



U-+(0) =- 



aT (1 + 6-/3) 

(1 - p + v^T^) 
, r(i-s-£) 
' r (1 - p - JWTtf) 



(57b) 



T tr (0.1) 



log(3). 



(61) 



Finally, although Eqs. (|5"9")) were obtained within the 
weak dephasing approximation, they are valid even in 
the strong dephasing regime. In this limit dynamics are 
completely incoherent and arc governed by rate equations 

M. 
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B. General case with jij — 7 

Although analytic solutions for Eqs. (f28|) cannot be 
derived when <fi(t) is no longer constant, the main fea- 
tures for the system dynamics can be qualitatively dis- 
cussed. In order to simplify Eqs. (|28|) . we first note that 
for Gaussian pulses the coupling term <f> sin 9 scales with 
the relative pulse delay r, see Eqs. (|65|) and (|66|). On 
the other hand the effective coupling terms f2 su (i), il sv (t) 
and Q uv (t) all scale with the relaxation rate 7. Then, in 
the weak dcphasing regime these latter coupling terms 
arc negligible compared to cb sin 6. Using this approx- 
imation the equations for the coherences u(t) and v(t) 
become 

u(t) « -T u (t)u(t) + 2v(t)<j)sin6, 

(62) 

v(t) w -T v (t)v(t) - 2u(t)0sin6>. 
Furthermore, we have that 

|r 5 .(t)|» \ci su [t)\,\ci sv (t)\. (63) 

With this the evolution of s(t) and consequently that of 
the dark states populations, is dominated by the effective 
relaxation rate r s (t), 

a(t) « -T,(t)a(t). (64) 

Thus, in the weak dcphasing regime coherences obey 
the dynamics of a two- level system Eq. (|6"2"|) . while the 
populations decay exponentially at a rate T s (t) (|64|) . 
Although analytic solutions could be derived for Eq. 
(f6"2"| following a similar method to that of Sec. IIV Al 
the complexity of the final expressions limits their use. 
Nevertheless, useful conclusions can be drawn by sim- 
ple inspection of these equations. When the relaxation 
rate 7 is increasing, and for fixed delays, the coherences 
and the populations will decay. Taking into account 
that r s (±oo) = r u (±oo) = 0, r„(±oo) = 7 and that 
fijj(±oo) = 0, wc anticipate that in the long time limit, 
v(t — > 00) = 0, whereas both s(t) and u(t) acquire a finite 
constant value. Thus in the long time limit coherences 
are partially preserved. 

The effect for different delay times on the system dy- 
namics is more complicated. This is due to the depen- 
dence that both the coupling <psm(6) and the effective re- 
laxation rates have with respect to the delay time. Tak- 
ing into account the results from the previous section 
IIV Al and those for STIRAP in A-systems we expect 
that the dependence of the fidelity with respect to the 
pulse delay will reflect the dependence of the transition 
time T tr (e) on the pulse delay. 

V. NUMERICAL SIMULATIONS 

We present now results from numerical simulations 
with the master equation (|16p and for Gaussian pulses. 



For overlapping control and Stokes pulses, f2 c (t) = 
we use the pulses given in Eq. (|44[) . For the pulse order- 
ing where the Stokes pulse proceeds and ends after the 
pump pulse, while the control pulse is delayed relative to 
both pulses, the parametrization reads 

n p (t) = noe-c*-/ 2 ) 2 /**, 

n s (t) = n e~ {t - T/2)2/{2T2 \ (65) 

ft c (t) = fVr (t+T/2)2/(2T2) . 

For the pulse ordering Stokes-control-pump we have 

fyt) = fV- (t+r/2)2/T2 , 

n s (t) = n a e- {t - T ^ 2 ' T \ (66) 

n c (t) = n e~ t2/T2 . 

For the ordering control-Stokes-pump, the pulses are the 
same as above with the only change that the Stokes and 
control pulses are interchanged, i.e. O s (t) — > f2 c (i) and 

n c (t) -> n,(t). 



A. Overlapping Stokes and control pulses 

Starting with a pulse ordering where the control and 
Stokes pulse overlap, i.e. f2 s (t) = fi c (i), in Fig. [3]we plot 
the populations pjjipo) for the bare states ipj, and for 
different relaxation rates 7^ = 7. The validity of the 
analytic results of Sc dIV Al is confirmed, where we see 
that the results almost coincide. The expected popula- 
tion damping for the dark states due to losses towards 
the bright states, is evidenced as a population drop for 
the states 1^3 and ^4. This is accompanied from popula- 
tion rise for the other two states 4>i and ^2 ■ We also note 
the predicted asymptotic behavior for strong dephasing, 
where pjj(oo) — > 1/4. 

Population losses for the dark states will also result in 
a fidelity decay for the final coherent state Eq. (jT5|) . For 
weak dephasing 7T <C 1 and for times shorter than I/7, 
i.e. T < t max < I/7, the fidelity ([30]) reads 

F 2 {t max ) = pl^oo) + Re{pl 2 (t max )}, (67) 

where Pii(oo) is given by Eq. (|59b|) and Re{p" 2 (*max)} 
is derived from Eq. (|59a|) . with the substitution 

lim e~ 7 ' — > e ~ 7tmax . 

t— >oo 

In Fig. 21 the fidelity is plotted against the relaxation rate 
-fij = 7 for the Gaussian pulses (|44|) . The dotted line is 
the fidelity F(t max ) for t max = 5T, where an exponential 
drop for increasing 7 is noted. More specifically, it can 
be shown that for 7T 2 /r <C 1 the fidelity is 

F 2 (t max )^ - (l + e-^j +± e -^-*»>: (68) 
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7 (units of T~') 

FIG. 3. (Color online) Final populations pjj(oo) for the bare 
states tpj plotted against the dephasing rate 7y = 7 for Gaus- 
sian pulses Eq. (|44|l . The peak Rabi frequency is Qo = 50T~ 1 , 
and the delay is r = 1.5T. The dots are numeric results and 
the solid lines were obtained from Eqs. (|27|l and (|59b[) . and 
with the aid of Eqs. fl8]l. (fi9)| and J2IJ. 
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FIG. 4. (Color online) The fidelity for the coherent superpo- 
sition state (|15|) as a function of the dephasing rate 7y = 7 
for the Gaussian pulses l|44p . The peak Rabi frequency is 
S7o = 50T~ and the delay is r = 1.5T. The dots are numeric 
results, whereas the solid line is the analytic result for strong 
dephasing or for t I/7. The dotted line is the fidelity for 
the weak dephasing regime (|67[) . and for t max — 5T. 



On the other hand for strong dephasing 7T 2 /r ^> 1, or 
in the long time limit t ^s> 1/7, the system completely 
decoheres and the only contribution to the fidelity is from 
the population ^(00) Eq. (|59bp . i.e. 

F 2 (oc)=PhM- (69) 

For this limit the final fidelity is well below unity, see 
solid line, and it rapidly drops for increasing dephasing 
rate. 

As already noted in Sec. IIV Al the delay r is expected 
to have an inverse effect on the fidelity. This is due to 
the inverse dependence of the transition time Eq. (|61[) 
with respect to the delay. This is demonstrated in Fig. [5] 
The fidelity for different peak Rabi frequencies is plotted 
against the delay for 7^ =7 = 0, Fig. [SJa), and for 
lij = 7 = T^ 1 , Fig. \5\b). From this we see that for 
adiabatic evolution the fidelity increases as Eqs. (|69|) 
and (|59b[) predict (dashed). The regime of validity for 
these equations increases for increasing pulse areas, and 
this is because of the adiabatic condition. 

B. Stokes-control-pump pulse ordering 

When considering a pulse ordering where the Stokes 
proceeds the control and the pump is delayed further 
from the control, Eq. (|66[) . the exact form of the result- 
ing coherent state (fT5|) will be a function of the delay r via 
the geometric phase i9 s ©. This means that the inverse 
effect that the increasing delay r has on the fidelity will 
reflect differently upon different coherent states. This is 
shown in Fig. [51 where the fidelity for different coherent 
states (different r) as obtained from numerical simula- 
tions is plotted against the dephasing rate 7^ =7. It is 



clear that as we increase the delay time and for a given 
7, the fidelity for the corresponding state 

*(oo) = -^sin^r) - ^ 4 cosier), (70) 

will be higher than the fidelity for a state that corre- 
sponds to smaller r. This, as already said, is because of 
the inverse dependence of the transition time pip on the 
delay r, see Fig. [7] As the relaxation rate increases, the 
fidelity for all the coherent states decreases reaching the 
strong dephasing asymptotic F = 0.5. 

C. Stokes proceeds and ends after the pump pulse 

In contrast to the previous two pulse orderings, this 
one differs in the effect that the delay r has on the fi- 
delity. Dephasing has the same effect where the fidelity 
drops exponentially for increasing 7^ = 7, see Fig. [8] In 
addition to this increasing the delay will also result in a 
fidelity reduction. The explanation for this is again given 
in terms of the transition time and its dependence with 
respect to the delay. 

As already pointed out in Sec. MI CI the fidelity for 
the coherent state (T!?)) and for t —> —00 is F 2 (— 00) = 
cos 2 i? £ ,. As the geometric phase increases, i.e. the de- 
lay increases, the initial value for the fidelity decreases. 
Thus, for increasing r the transition time is expected to 
increase too. This is shown in Fig. |H1 where the fidelity as 
a function of time and for different delays in the absence 
of decohcrence, Fig. [H^a), and the corresponding transi- 
tion times, Fig. |UJb), arc plotted. As we can see in Fig. 
[5{a) the initial value for the fidelity reduces for increasing 
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FIG. 5. (a) The fidelity for the coherent state (fT5|) plotted 
against the delay time r for the Gaussian pulses (53) , as 
obtained from simulations with the master equation. The 
peak Rabi frequencies fi are denoted next to the respective 
curves and jij =7 = 0. (b) The same with (a) but for 
7ij = 7 = T -1 . The dashed line is the analytic result ob- 
tained from Eqs. ((69j) and (|59b|) . 



FIG. 6. (Color online) The fidelity for the coherent state (|13[) 
as a function of the dephasing rate 7^ = 7 for the Gaussian 
pulses (|66l) . The peak Rabi frequency is flo = 200T -1 and 
the delay is r = LOT (solid) r = 1.5T (dashed) and r = 
2.0T (dotted). The results were obtained from a numerical 
simulation with the master equation Q5) . 




Delay (units of T) 



delays. For this reason the transition time will increase, 
Fig. [UJb), and consequently the effect of dephasing on 
the fidelity increases. 



FIG. 7. The transition time Tt r (e) (|31[1 for the Gaussian 
pulses (|66|l . The peak Rabi frequency is fio = 200T -1 and 
e — 0.1. The times £ e and ii_ e , were derived by numerically 
solving the equations F 2 (t e ) — e and F 2 (ti_ e ) = 1 — e. 



VI. CONCLUSION 

In this paper we have studied the effects that dephas- 
ing has on STIRAP in tripod configurations. We have 
derived an exact adiabatic solution for the master equa- 
tion in the case of overlapping Stokes and control pulses, 
and for weak dephasing. The results were verified with 
numerical simulations to provide a very accurate approx- 
imation for the fidelity dynamics. In the adiabatic limit, 
population losses and dephasing for the dark states de- 
pend only on the relative relaxation rates of the three 



bare ground states, but not on those for the intermediate 
state. 

The fidelity exponentially decreases for increasing de- 
phasing, whereas the pulse delay has an inverse effect. 
This is due to the fact that the transition time decreases 
for increasing delay. This way dephasing effects are sup- 
pressed. Using numerical simulations we extended our 
studies to different pulse ordcrings. For a pulse ordering 
of the form Stokes-control-pump, or when interchanging 
the order of the control and Stokes pulses, similar dynam- 
ics are observed as with overlapping Stokes and control 
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FIG. 8. (Color online) The fidelity for the coherent state (JT2J) 
as a function of the dephasing rate <y<j = 7 for the Gaussian 
pulses (|65l) . The peak Rabi frequency is fio = 200T" 1 and 
the delay is r = 0.5T (solid) r = LOT (dashed) and r = 
1.5T (dotted). The results were obtained from a numerical 
simulation with the master equation (| lCjf) . 



pulses. The fidelity decrease exponentially with the de- 
phasing, whereas the delay has again an inverse effect. 
This will reflect differently on each coherent state. This 
is because the geometric phase, which characterizes the 
final superposition, is a function of the delay. 

For a pulse ordering where the Stokes pulse proceeds 
and ends after the pump pulse, while the control pulse 
follows the delay has a different effect. The reason for 
this is that the transition time increases with the delay. 
Dephasing has again the same effect leading to an expo- 
nential decrease for the fidelity. 
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FIG. 9. (Color online) (a) The fidelity for the coherent state 
(|12[) as a function of time for the Gaussian pulses (|65[) . as 
obtained from simulations with the master equation. The 
peak Rabi frequency is fio = 200T -1 , whereas the delay is 
t = 0.5T (solid) r = LOT (dashed) and r = 1.5T (dotted). 
The relaxation rates are 7;-,- =7 = 0. (b) The transition 
time Tt r (e) Eq. (32} for e = 0.1. The times t 1+€ and ti- e , 
were derived by numerically solving the equations F 2 (ti+ € ) = 
(1 + e) cos 2 (tf 9 ) and F 2 (ti_ e ) = 1 - e. 



Appendix A: Deriving the equations for the 
two-level system 

From Eqs. CTZD, flU]), (JTSJ), flU, and (25) we can derive 
the following equation for the population inversion W34 = 

Ph ~ P44 

w 34 (t) = (712 sm 2 {9) + cos 2 (0)r 2 (0)) w 34 (t), (Ala) 
where 

Tj ((/)) = cos 2 (0) 7j3 + sin 2 (0) 7j - 4 . (Alb) 

Taking into account the initial condition p 33 (— 00) = 
p2i(— 00) = 0, or W34(— 00) = 0, we see that that the 



inversion is w 34 {t) = 0, i.e. P33 = /o| 4 . Making next the 
substitution 



P33 = P44 = (1 - Pii P22)A 



(A2) 



we get the following master equation for the effective two- 
level system spanned by the two dark states @ 

Ph = -V%f% - ij>sm(e)[a z ,p d ] kl - D° kl , (A3a) 
where p d is 

^ -{%%)• (A3b) 
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and cr 2 is the relevant Pauli matrix. In this equa- 
tion the Einstein summation convention is used, where 
i,j,k,l = 1,2. The tensor X>jy and the matrix re- 
sulted when using the rotation matrix R(t) (|19p to trans- 
form the master equation (|16[) in the adiabatic basis Eq. 
(|20[) . At this point their exact analytic form is not impor- 
tant. Instead some very useful properties that are going 
to be used next are listed below 



and 



T) 11 — T) 22 
u 2\ — u ' 



12' 

T) 11 — T) 22 

21: U \2 — U VXl 



7? 11 — T) 22 

u \\ ~ u 22i 



jy22 _ jyll 



■t"\12 T")12 T")21 T")21 

u \\ ~ ^22; u \\ ~ ^22 > 

v\l = {v\ly = v^ = {vl\)\ 

Re{D{ 2 } = ^Re{VH}, 
D°n=D° 22 = -^1+^1), 

D°i2 = (D° 21 y = -\v\\. 



(A4a) 
(A4b) 
(A4c) 
(A4d) 
(A4e) 

(A4f) 
(A4g) 
(A4h) 
(A4i) 



Having these properties we can now proceed and de- 
rive the equations for the population inversion w1 2 {t) = 
Pii ~P22i the populations p\ x and p 22 , and the coherences 
u(t) = V2Re{pl 2 } and v{t) = \/2Im{pJ 2 }. Starting from 
the inversion wf 2 (t), is easy to show that 

< 2 (t) = (»n-»uW 2 (i)- (A5) 
Taking into account the fact that 6(— oo) = and 
Pli(— oo) = 1, we have that p"-^— oo) = P22^ 00 ) = 
and iof 2 (- oo) = 0. From the above equation we have 
that at all times 



«>£,(*) = o 



and that 



P22 - 4 



2 ' 



(A6) 



(A7) 



With this result the remaining three equations for s(t), 
u(t) and v(t) (|28|) are derived, with the initial conditions 
being s(— 00) = —1/2, u(— 00) = l/\/2 and v(— 00) = 0. 
The effective relaxation rates and Rabi frequencies read 



r s (t) = p 1 1 1 1 + p22 

= 1 sin 2 (20)r 1 (0) + i cos 4 (#) sin 2 (20)734, (A8a) 
= i sin 2 (20)r 1 (0) + i (1 + sin 2 (0)) 2 sin 2 (20)734, 



r„(t) = 2>J§ - Re{X> 2 12 } 

= cos 2 (6>)r 1 (0) +sin 2 (6>) cos 2 (20)734, 



(A8b) 
(A8c) 



n«,(t) =Re{©^} - -sin 2 (20)r 1 (0) 

- - cos 2 ((9) (1 + sin 2 ((9)) sin 2 (20) 7 34, (A9a) 
Cl sv (t) =Im{V\l} = -icos 2 (0)sin(0)sin(40) 7 34 

+ - cos 2 (0) sin(0) sin(20)(7i4 - 713), (A9b) 
=Im{2?i 2 } = (sin(30) - 7sin(0)) sm(40) 73 4 

- - cos 2 (0) sm(6) sin(20)(7i 4 - 713). (A9c) 



Appendix B: Adiabatic integrals 

In order to obtain s(oo) and u(oo), we first use the 
inverse rotation 7£ _1 (oo) ([38]) on the state ip(oc) (|50|) 
to get ci(oo) and 02(00) (|36p . Using the expressions for 
c±(oo) and Eq. (f36|) we have for s(oc) 



s(oo) oc exp(/ s ) 
where is 



(Bla) 



/ (e + (t)-T s (t))dt+ / (e_(t)-r a (t))di 
J -oo Jo 

(Bib) 

and for m(oc) 

it(oo) oc exp(I ll ) = exp (J (e + (t) - T s (t))dt 

(Blc) 

Although the expressions for e±(t) and r s (i) are very 
complicated, they can all be parametrized in terms of 
single dimensionless variable x = Atr/T 2 , i.e. 



e±(£) = 79±(x), 
and 

T s (t) = 73s 0), 
where 

S-O) 



(B2a) 



(B2b) 



4(6^ - 1) 



(1 + e*)(2 + e*) 2 (l + y/l + 8/(1 + e x ) 2 ) ' 

(B2c) 



, s (l-e 2 -)(l + Vl + 8/(1 + e^) 
3+0) = 0/0 1 ^2 ' ( B2d ) 



2(2 + e x ) 
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and 



Breaking the integral into three parts we have 



9s(x) 



2(l + 2e x ) 
(2 + e*) 2 ' 



(B2e) 



With this the integral for I s in Eq. (|Bla[) . takes the 
form 



7T 2 f° 

I s =— (g+(x) -g s (x))dx 

J — OO 

7T 2 [°° 
+ J (g-(x) - g s {x))dx, 



(B3) 



where both integrals converge, and can be calculated nu- 
merically giving 



4r 



(B4) 



where c s — 2.42. When calculating the integral I u in Eq. 
(IBlbl) . attention must be paid for x — > 00. For this limit 
the integrand converges to a finite value 



7T 2 f° 

Iu=— {g+(x) - g s (x))dx 



ill 
4r 



2 roc 



(g+(x) - 9s{x) + l)dx 



(B6) 



ill 
4r 



2 r x ~ * 00 



dx, 



where in the integrand for x > we have added and 
subtracted its asymptotic value for x — > 00. The first 
two integrals converge, whereas the third one gives rise 
to the exponential term e ft in Eq. (|59a[) . Thus the 
integral I u reads 



CujT 2 
4r 



(B7) 



lim (g+{x) -g s (x)) = -1. 



(B5) where c u = 0.68. 
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